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Abstract 

We consider time dependent correlation functions of non-Abelian gauge fields at finite 
temperature. An effective theory for the soft (p ~ g 2 T) field modes is derived by in- 
tegrating out the field modes with momenta of order T and of order gT in a leading 
logarithmic approximation. In this effective theory the time evolution of the soft fields 
is determined by a local Langevin-type equation. As an application, the rate for hot 
electroweak baryon number violation is estimated as T ~ g 2 log(l/g)(g 2 T) 4 . Further- 
more, possible consequences for non-perturbative lattice computations of unequal time 
correlation functions are discussed. 
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Finite temperature field theory for non-Abelian gauge fields is non-perturbative for 
soft[| momenta p ~ g 2 T, where g is the gauge coupling and T is the temperature 
|l|, . The dynamics of the soft gauge fields determines the rate for electroweak baryon 
number violation at very high temperatures ||. This rate is determined by a real time 
correlation function of the type 

C(t 1 -t 2 ) = {0[A(t 1 )]0[A(t 2 )}) (1) 

where (• • •) denotes the thermal average and C?L4(£)] is a gauge invariant function of 
the gauge fields A a {t,-x). 

The aim of this letter is to derive an effective theory for the dynamics of the soft 
gauge fields by integrating out the "hard" (k ~ T) and "semi-hard" (k ~ gT) fields. 
For weak gauge coupling g this can be done perturbatively. This effective theory should 
allow for a non-perturbative computation of real time correlation functions like ([!]), e.g., 
on a lattice. 

For time-independent physical quantities, like the free energy or correlation lengths, 
the effective theory for soft gauge fields is well established. It is a three-dimensional 
gauge theory, the parameters of which are determined by dimensional reduction 0, 
||. The dimensionally reduced theory can serve as an input for Euclidean lattice 
simulations or other non-perturbative methods 0. The three-dimensional theory 
is much simpler to simulate on a lattice than the full four- dimensional one. 

For time dependent quantities it may be even more important to find an effective 
theory for the soft field modes: There is no apparent way to calculate real time quanti- 
ties like (0) in lattice simulations of finite temperature quantum field theory. Presently 



the only known tool to evaluate them is the classical field approximation ||-[10] and 



variants thereof |Ll|-||14| which contain additional degrees of freedom representing the 
hard field modes. 

The reason why the classical field approximation is expected to be reliable is that the 
field modes with soft momenta have a large occupation number and should therefore 
behave classically. The high momentum modes have occupation number of order unity 
and do not behave like classical fields. They are however weakly interacting and can 
be treated as almost free massless particles moving quasi-classically. 



There has been a long discussion (see Refs. |15 , 16fl ) about whether these hard par- 
ticles do affect the dynamics of the soft field modes, and, in particular, whether they 
play a role in electroweak baryon number violation. If all field modes with momenta 
larger than g 2 T were irrelevant, there would be only one scale in the problem and on 



2 For spatial vectors we use the notation k = |k|. 4-momentum vectors are denoted by if = (k°, k) 
and we use the metric K 2 = k 2 — k 2 . P always refers to a soft momentum and K to a semi-hard or a 
hard one. 
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Figure 1: Sub- leading contributions to the effective theory for the soft (p ~ g 2 ? 1 ) gauge 
fields. The loop momenta are semi-hard (A; ~ gT). The small dots are "bare" (non hard 
thermal loop) vertices. The full lines denote hard thermal loop resummed propagators, 
the dashed lines represent ghost propagators. 



dimensional grounds the rate can be estimated as F ~ {g 2 T) 4 ||16|| . Arnold et al. JT7] 



demonstrated that the dynamics of the soft fields is damped by the hard particles and 
they obtained the estimate V ~ g 2 (g 2 T) 4 . There are, however, contributions to the soft 
dynamics due to semi-hard field modes which are as important as the hard particles 
im|. These contributions are the subject of this letter. Here only the main points of 



the calculation will be described. Further details can be found in Ref. [18]. 

The first step in deriving an effective theory for soft gauge fields is to integrate out 
the hard field modes with spatial momenta of order T. The dominant contributions 



from these fields are the so called hard thermal loops |Tj|. A consistent perturbative 



expansion for the semi-hard field modes requires the use of the so called Braaten- 
Pisarski scheme, or hard thermal loop effective theory^. In this scheme the hard thermal 
loop propagators and vertices must be treated on the same footing as their tree-level 
counterparts. 

In the present context integrating out the hard field modes yields an effective theory 
for the semi-hard and the soft fields. As we will see below, the leading contributions 
to the effective soft dynamics arise from hard thermal loop induced interactions of the 
soft fields with the semi-hard ones. Since there are no hard thermal loop vertices in 
Abelian theories, the terms we are going to compute do not have an Abelian analogue 



(cf. Refs. 0, [22 



In order to integrate out the semi-hard fields we introduce a separation scale \i such 
that 

g 2 T < (i < gT (2) 

and integrate out all field modes with k > /i. 

Let us first discuss what kind of diagrams are relevant to our problem. We have to 

consider diagrams in the hard thermal loop effective theory in which all internal loop 

3 For certain momenta K — (fc ,k) with k of order gT, however, the hard thermal loop effective 
expansion breaks down, e.g., when K is on the light cone pi. 
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Figure 2: Leading contributions to the effective theory for the soft (p ~ g 2 T) gauge 
fields. The loop momenta are semi-hard (k ~ gT). The full lines denote hard thermal 
loop resummed propagators and the heavy dots are hard thermal loop vertices. 

momenta are semi-hard and the external momenta are soft. These have to be compared 
with the corresponding "tree- level" terms like, i.e., the transverse hard thermal loop 
polarization operator which is given by 

SU t (P) = ijrf(p) ml J ^V^. (3) 

Here Pt J {p) is the transverse projector 

if (P) = S* - % (4) 

and m\y is the leading order Debye mass squared. In a pure SU(iV) gauge theory, it 
is given by m 2 ^ = (l/3)Ng 2 T 2 . In the hot electroweak theory it receives additional 
contributions due to the Higgs and fermion fields. Furthermore, = (l,v), and the 
integral J dQ v is over the directions of the unit vector v, |v| = 1 . 

Consider the corrections to Eq. (|3|) due to the semi-hard field modes. The leading or- 
der behavior of the diagrams in Fig. [I] can be easily estimated since they are completely 
analogous to hard thermal loop diagrams for semi-hard external momenta except that 
now all momenta are smaller by a factor g. The k integral behaves like / dk ~ gT 
which gives a suppression factor g relative to the hard thermal loops. In other words, 
the diagrams in Fig. [I] can be neglected relative to the "tree- level" term (^). 

Now consider the analogous diagrams (Fig. |2|) with the "bare" vertices replaced by 
hard thermal loop vertices. There are only two such diagrams since there is no hard 
thermal loop vertex involving ghost fields. The vertex in Fig. 0(a) contains a term 
proportional to g 2 po /(v-P) 2 times a factor of order gT. Here the four vector is of the 
same type as in Eq. (|3]). Due to the Bose distribution function in the loop integral there 
is a factor T so that diagram 2(a) can be estimated as m 2 :) g 2 Tpo/(v -P) 2 . Comparing this 
expression with Eq. (^) we see that both terms are of the same order of magnitude. In 
fact, the diagrams in Fig. |2] can be even larger by a factor of \og(gT / //): The transverse 
propagator is unscreened when |fc | *C k and in the infrared the loop integrals can be 
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Figure 3: Contribution to the effective theory for the soft (p ~ Q 2 T) gauge fields 
involving hard thermal loop n-point vertices. The notation is the same as in Fig. |2|. 

sensitive to the separation scale \x. We will see below that this logarithm indeed occurs. 
It is specific to the transverse semi-hard field modes. In this letter we will compute 
only these leading logarithmic contributions to the effective soft dynamics. 

In order to obtain the effective theory for the soft field modes we have to consider 
diagrams with more external soft fields and with more internal semi-hard propagators 
like in Fig. |3|. Higher hard thermal loop n-point vertices contain a large number of 
terms which make a diagrammatic analysis very difficult. It is more convenient to 
use an effective field theory description and solve the corresponding field equations of 
motion. 

There are locaf] formulations of the hard thermal loop effective theory due to Blaizot 
and Iancu |2^| and due to Nair |[24|| . These formulations make the physical content of 
the hard thermal loop effective theory quite transparent (see also Ref. fl25|). The 



formulation due to Blaizot and Iancu is the non-Abelian generalization of the Vlasov 
equations for a QED plasma. It consists of Maxwell's equation for the semi-hard and 
soft gauge fields, i.e., fields with spatial momenta of order gT or less, 

[D„,F^(x)]=f(x) (5) 

where = — igA^x) is the covariant derivative and A^{x) = A a ^(x)T a with 
hermitian generators normalized to Tr(T a T b ) = (l/2)S ab . The current j u (x) is due to 
the hard field modes. These modes are weakly interacting and they behave as massless 
particles moving at the speed of light with 3-velocity v. The current can be written as 

f{x)=ml j^v-W^v). (6) 

The field W(x, v) = W a (x,v)T a is proportional to the deviation of the distribution of 
hard particles from the equilibrium distribution. It is determined by the equation of 
motion 

[v-D,W{x,v)] = v-E(x), (7) 



*The hard thermal loop vertices are non-local in space and time. 
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where E is the electric field strength. The conserved Hamiltonian corresponding to 
Eqs. (D-(0) is @ 

H = J d 3 xTr |e(z) • E(x) + B(x) • B(ar) + ml J ^W(x, v)W(x, v) j . (8) 
Originally, the hard thermal loop effective theory was designed for calculating Green 



functions with semi-hard momenta. Following |1 1] we assume that the leading order 
dynamics for both soft and semi-hard fields is correctly described by the hard thermal 
loop effective theory. More precisely, we assume that real time correlation functions of 
gauge invariant operators like (HD can be obtained as follows: 

• Compute the solution A a Jx) of the non-Abelian Vlasov equations ©-© for given 
initial conditions v4 in (x), E in (x) and Wi n (x, v) at t — 0. 

• Then C(ti— t 2 ) is given by the product 0[y4(ti)](9[y4(t 2 )] averaged over the initial 
conditions with the Boltzmann weight exp(— (3H[A- m , E; n , Wi n ]) where Gauss' law 
has to be imposed as a constraint. 

In the average over initial conditions ultraviolet divergences occur. These are due to 
the fact that a classical field theory at finite temperature is not well defined. One has 
to use an UV cutoff A which satisfies 

gT < A < T. (9) 

After the semi-hard field modes have been integrated out, the dependence on A cancels 
against the appropriate A-dependent counter-terms to be included in the hard thermal 
loop effective theory^. 

To integrate out the semi-hard modes, we decompose A, E and W into 

A -> A + a 
E -> E + e 

W -> W + w. (10) 

The soft modesQ A, E and W contain the spatial Fourier components with p < \x while 
the semi-hard modes a, e and w consist of those with k > p. Let us note that the 



separation flip ) does not respect gauge invariance. Nevertheless, the results we are 



going to obtain, will not depend on the choice of a gauge. Furthermore, we do not 



5 See also the discussion in Ref. Q 

6 For notational simplicity we do not introduce new symbols for the soft modes. From now on A, 
E and W will always refer to the soft fields only. 
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specify the way the cutoff is realized precisely. The only cutoff dependence which we 
will find is a logarithmic one for which the precise choice should be irrelevant. Here we 
consider only the transverse semi-hard gauge fields. 

Due to the non-linear terms in Eqs. (||) and (^) the equations of motion for the 
soft and the semi-hard fields are coupled. Let us see which of these couplings are 
relevant to our problem. The non-linear terms in Maxwell's equations (|5|) correspond 
to non-hard thermal loop vertices. For the interaction between soft and semi-hard fields 
these vertices can be neglected as we have already argued. The important vertices are 
the hard thermal loop ones. They correspond to non-linear terms in Eq. (]?]). Thus, 
substituting Eq. (|10D into Eqs. @-(Q), the equations for the soft fields read 

[D„F^(x)]=ml J ^±L v »W(x,v) (11) 

and 

[v ■ D, W{x, v)} = v • E(z) + f (ar, v), (12) 

where 

e(*,v) =/ c (vaJ(xK(^)) sof . (13) 

Here the subscript "soft" indicates that only spatial Fourier components with p < jj, 
are included. The structure constants f abc are defined such that [T a ,T b ] = if abc T c . 
The transverse gauge field a t (x) and the field w(x,v) in Eq. fll3|) are to be understood 
as the solutions to the coupled equations of motion (pT|)-(p^) and (|14l)-(|16|). 

In the following we discuss how the semi-hard fields can be eliminated from Eqs. 
(|TT|)- ([T3|) leading to the final result (p9|). A more detailed description of the calculation 



can be found in Ref. [flSfl . 

The semi-hard fields are weakly interacting and can therefore be treated perturba- 
tively. Since we are only interested in leading order results, it might appear sufficient 
to insert the solutions of the free equations of motion, a t0 (x) and w (x, v) into £(x, v). 
The corresponding term will be called £,o(x, v). As will be seen below, a next-to-leading 
order term has to be included as well. For this term, which will be denoted £\(x, v), 
one can neglect those interactions in Eq. (0) which contain only semi-hard fields. In 
the diagrammatic language such interactions correspond to parts of hard thermal loop 
vertices which do not contain the enhancing factors 1/ (v ■ P) discussed above. The rel- 
evant non-linear terms are the ones which contain a product of a soft and a semi-hard 
field. Thus the equations of motion in spatial momentum space can be written as 

a a t (t,k) + k 2 a a t (t,k) = mlj^v t w a (t,k,v), (14) 
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w a (t, k, v) + zv • kw a (t, k, u) = v • e t (t, k) + /i a (t, k, u), (15) 

where v\ = P^ J (k)v J is the transverse projection of v. The term h a (t, k, t>) contains the 
interaction of the semi-hard fields with the soft ones. In coordinate space it reads 

h a (x, v) = gf abc (-v ■ A b {x)w c {x, v) + v • a%(x)W e (x, vj) . (16) 

A convenient method for solving Eqs. (]T^)-([TBD is the Laplace transformation^. The 
Laplace transform of a function f(t) is defined as 

roc n 

f(k°) = / dte ik0t f(t) (17) 
Jo 

and it is an analytic function in the upper half of the complex fc°-plane. Applying (|TT| ) 
to Eqs. ( |HD and ( Jl5| ) one obtains 

al a (K) = a™(K) + J ^A\ 2 (K,v 1 )h a (K,v 1 ) (18) 

and 



w a (K,v) = w a (K,v)+ I ^Al 2 (K,v,v 1 )h a (K,v 1 ), (19) 

where & t o{K) and wq{K, v) are the solutions to the free equations of motion. They are 
determined by the initial values at t = 0. The propagator functions in (|Tg|), flHf ) are 
given by 

A\ 2 (K,v) = iml^G^K)^—, (20) 

v ■ K 

A 22 (K,v, Vl ) = W s2 )(v- Vl )^--m^XGt J (^) f 77- (21) 

Here G^(JC) is the transverse hard thermal loop resummed propagator, 

GV(K) = 5-^ — Pl j (k). (22) 

* v ; -K 2 + m t (K) K J K J 

Furthermore, f^ 5 ^ is the delta function on the two dimensional unit sphere: 

/^ / (v 1 )^(v-v 1 )=/(v). (23) 

The solutions to the equations of motion, Eqs. ( |T8"D and (|H|), are still formal in the 
sense that h a (K, v) contains the complete solutions a t (iT) and w(K, v) as well as A(P) 
and W(P, v). However, by iterating Eqs. ( |TB[ ) and flT5| ) one obtains a series in which 

7 Or one sided Fourier transformation, see,e.g., Ref. pjl. 



each term contains only the free solutions ato(K) and wo(K,v) together with the full 
A(P) and W(P,v). Since Eqs. (|I4D and fllBD are linear in a t and w, the solutions are 
linear in a t0 and wq- The function £(x,v) is therefore bilinear in the fields a t0 and Wq 
which always appear in the form 

X = <f>(K-P)(f>'(P'-K), (24) 

where the <ft and <$' are either a t o or wq. The solution to the (non-perturbative) equa- 
tions of motion for the soft fields will contain products of x' s - To obtain a correlation 
function like ([I]) one has to average over the initial conditions for both soft and semi- 
hard fields. 

In general, the thermal average of a product Xi' "Xn can be approximated by discon- 
nected parts: The connected part (xi • • • Xn)c has one momentum conserving delta func- 
tion for a soft momentum. A disconnected part like for instance (xi • • ■ Xm)c(Xm+i 4 4 4 Xn)i 
contains two such delta functions. The corresponding contribution to a non-perturbative 
correlation function has one soft 4-momentum integration less and one additional semi- 
hard 4-momentum integration. We can therefore estimate^ 

(Xl 4 • • Xn) c ~ 9 S -^(XI 4 4 4 Xm) c (Xm+l 4 4 4 Xn) ^ (25) 

Here p° is the typical frequency scale for the soft non-perturbative dynamics. This 
scale will be estimated from the equations of motion for the soft fields which we are 
going to derive (see Eq. fl45f) ). This consideration suggests that £(x,v) in Eq. ( |12D 



can 



be simplified by replacing the x% by their thermal averages: x% ~^ (Xi)- For the lowest 
order term £ , which does not depend on the soft fields, this replacement gives zero due 
to the antisymmetry of the structure constants f abc . Therefore we have to leave £o i n 
Eq. (|TJ) as it stands and the solution to the equations of motion for the soft fields will 
contain products of £q's. The thermal averages of these products can be approximated 
by a product of two point correlators. We also have to include the next-to-leading 
order term £1 which is linear in the soft fields A and W. 

Neglecting the external soft 4-momentum inside loop integrals the leading logarith- 
mic result for the two point correlator of £o i n configuration space reads 

%(x 1 ,v 1 )e (x 2 ,v 2 )) = -2iV^f log ( 9 —\ I(v 1 ,v 2 )5 ab 5 {4 \x l -x 2 ), (26) 



with 



f(,,»,)E-iW(T-»j4 / (v ; Vi) . (27) 

1 — (v • vi r 



7l 2 



8 As we have already discussed, the integration over semi-hard momenta can give logarithms of the 
separation scale [i. These logarithms are not included in this estimate. 
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The logarithmic /i- dependence in Eq. (p6|) arises because the transverse propagator 
(p!|) is unscreened for |A; | <C fc. Non-logarithmic contributions have been neglected 
in Eq. (|26|). To compute them, one would have to take into account the longitudinal 
semi-hard gauge fields as well. Since the external soft momentum has been neglected 
in order to obtain Eq. (p6|), this expression is valid only when applied to problems with 
a relevant time scale much larger than (gT)^ 1 . 

For £1 we can replace the products of the free semi-hard fields a t o and wo by their 
expectation values. Then the terms containing A(x) cancel. With the same approxi- 
mations which were used for Eq. (|26|) one obtains 



Ct(x,v) = Ng 2 T\og J ^Ll( v , Vl )W(x, Vl ). (2* 



Note that without the cancelation of the A(x)-dependent terms Eq. ( p9|) would not be 
gauge covariant. 

With the above approximations the equation of motion fli~2l ) reads 



[v ■ D, W(x, v)] = v • E(x) + £ (x, v) + Ng 2 T log (^j J 



{ -I(v,v 1 )W(x,v 1 ). 

(29) 



Let us see whether Eq. ([29]) is consistent with Maxwell's equation (|TTJ) for the soft 
fields which requires the current on the r.h.s. of Eq. (|TT| ) to be covariantly conserved. 
Integrating Eq. (|29| ) over the direction of v the third term on the r.h.s. drops out due 
to 



47T 



■l(v,v 1 ) = (30) 



and one obtains 



[D„f{x)]=mi I ^-Mx,v). (31) 



That is, the current appears not to be conserved. However, as we have argued, only 
the two point function of £o should be relevant for the leading order behavior of a soft 
correlator. But the two point function of £0(^15^1) with the r.h.s. of Eq. fl3T|) vanishes 
due to Eq. (|30|). Thus we can replace 

-Ux,v)^0 (32) 



47T 

so that the current is covariantly conserved within the present approximation. 

Let us further note that, due to the term £ 0) Eq. ( |2^ ) appears to be non-covariant 
under gauge transformations. Nevertheless, a gauge invariant correlation function com- 
puted via Eqs. fl2HD and Q2DD will not depend on the choice of the gauge after the thermal 
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average over initial conditions since the two point function (^6|) is invariant under gauge 
transformations of £o- 

Eq. fl29|) is a Boltzmann equation for the soft fluctuations of the particle distribution 
W(x,v). The r.h.s. contains a collision term which is due to the interactions with the 
semi-hard fields. When diW on the l.h.s. of Eq. ( f29|) is of order (g 2 T)W, the collision 
term is larger than diW by a factor \og(gT//i). The reason is that the mean free path 
for the hard particles interacting with fields with k > fi is of order (g 2 T log(gT / 
which is smaller than the soft length scale (g 2 T) _1 by the same logarithm. The collision 
term is accompanied by the noise term £ which is due to the thermal fluctuation of 
initial conditions^ of the fields with k > fi. 

For a QED plasma there is no collision term at this order in the coupling constant. 
In this case the size of the collision term is determined by the transport cross section 
which corresponds to a mean free path of order order (e 4 T) _1 (cf. the discussion in 
Ref. For a non-Abelian plasma the relevant mean free path is determined by 

the total cross section which is dominated by small angle scattering: Even a scattering 
process which hardly changes the momentum of a hard particle can change its color 
charge which is what is seen by the soft gauge fields. 

What can one learn from Eq. (|29|) about the time scale relevant to non-perturbative 
physics? Let us first see how many field modes can be integrated out perturbatively. 
The mean free path for the hard particles interacting with the modes with k > fi 
is of order (g 2 T log(gT//z)) _1 . Perturbation theory for the hard particles must break 
down at this length scale. Thus we can decrease \i until \i ~ g 2 T ' log(gT / /i) . Then the 
logarithm in Eq. (f29|), in a first approximation, becomes log(l/g). 

We will now simplify Eq. (^9|) by neglecting terms which are suppressed by inverse 
powers of log(l/g). We introduce moments of W(x,v) and £, (x,v), 

W(x) = J^L W (x,v), (33) 

W h - in (x) = l^v h ---v in W(x,v), (34) 

""•(•'•) = J^v^--.vHo(x,v). (35) 

Multiplying Eq. (|29|) with the appropriate factors v 11 ■ ■ ■ v %n and integrating over the 
direction of v one obtains a set of coupled equations for the moments of W(x, v). From 
now on we will use the temporal axial gauge A = which is the most convenient for 
our purpose. The zeroth moment of Eq. ( p9|) is the equation for current conservation 



9 The term (o(i,w) should not be confused with the stochastic force discussed in Refs. [jlT], ^7], |30| |. 
The latter is due to the fluctuations of the initial conditions for the soft W(x, v) and it is also present 
in the theory described by Eq.(E9h. 
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(j32j). Taking the first moment of Eq. (|29"D gives 

d W%x) - \D\ W*{x)\ = ~E%x) + H{x) - \og(l/g)W i (x). (36) 

The l.h.s. of Eq. ( |36| ) is logarithmically suppressed^ relative to the term oc W l (x) on 
the r.h.s. and can be neglected. In this approximation we can determine W t (x) in 
terms of E(x) and £,q(x) without solving any differential equation. Introducing 

and 

4tt ml 1 

7 3 Ng*T\og(l/gy { ) 

the spatial components of Eq. (|ll|) become 

- E i + [Dj, F ji (x)} = 1 E i (x) + C{x). (39) 

The component, or Gauss' law, now reads[^] 

[D i ,E i (x)]=j (x). (40) 

The charge density jo(x) = m 2 :) W(x) satisfies 

d ja(x) = -yj (x)-[D i ,t i (x)]. (41) 



Eqs. (|3"9"D-(f4"l"D form a closed set of equations for the gauge fields and the charge density 
jo(x). The stochastic force C(x) (cf. Eq. (^)) satisfies 

(( ia (xi)( jb (x 2 )) = 2T 1 5 ij 6 ab 5^(x 1 - x 2 ). (42) 



Eqs. (p9|)-(^2|) are gauge covariant Langevin-type equations. The dynamics of the soft 
fields is Landau-damped by the hard field modes (see, e.g., Ref. pSfl ). The kinematics 



of the hard modes being massless particles moving on straight lines no longer appears in 
these equations because the semi-hard fields randomize^ the color of the hard particles 
on a length scale (log(l / g)g 2 T)~ 1 which is small compared to the length scale (g 2 T)~ 1 
of non-perturbative physics. The random force ( m ( x ) is due to the thermal fluctuations 
of the initial conditions for the semi-hard gauge field modes and for w(x, v). 
10 Provided that d <g 2 T, see below. 

11 For the charge density in Eqs. (|4C|)-([Ilj) we use the same symbol as for the one in Eqs. (|5|)-(^|).The 
latter has both soft and semi-hard Fourier components. 
12 I thank Guy Moore for this explanation. 
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We will now solve the linearized form of the equation of motion (p9|) . Even though 
we are interested in the non-perturbative dynamics of the soft gauge fields, this will 
allow us to determine the frequency scale on which the dynamics of these fields becomes 
non-perturbative. The solution for the transverse field reads 

A*(P) = -- P 2 [G t (P) - G t (0,p)]^ in (p) + G t (P) f-^(p) + CW] , (43) 
Po 

where the subscript "in" refers to the the initial values at t = 0. The transverse 
propagator G t is given by 

G t (P) = p2 l . ■ (44) 

-P 2 - I'fpo 

The typical size of the soft field is A t (:r) ~ A t> i n (x) ~ gT. Thus both terms in a 
covariant derivative <9 M — igA^ are of the same size making the thermodynamics of the 
soft fields non-perturbative. 

The dynamics of the soft gauge field becomes non-perturbative if it changes in time 
by an amount AA which is of the same size as A itself. Then the non-linear terms in 
the equations of motion are as important as the linear ones. AA t can be estimated 
as AA t (p) ~ p 2 G t (P)A tiin (p). From Eq. (g|) one reads off that AA t (p) ~ A t , in (p) 
when po ~ g A log(l/ 'g)T '. Thus we find that large non-perturbative changes AA t of the 
transverse (magnetic) gauge fields are associated with the frequency scale 

po~0 4 log(l/0)r. (45) 

Let us now use this result to estimate the rate for electroweak baryon number vio- 
lation at very high temperatures ||, |0J-[]T7|], i.e., well above the critical temperature 



T c ~ 100 GeV of the electroweak phase transition or crossover. At such high tempera- 
tures the Higgs field acquires a large thermal mass and decouples from the dynamics. 
Then it is sufficient to consider a pure SU(2) gauge theory. Baryon number noncon- 
serving processes are due to topology changing transitions in the electroweak theory 
which are characterized by a change AiVcs = ±1 of the Chern-Simons number 

Ncs = -£^e ijk J d*x (f«AI - 9 -e abc A*A)Al} . (46) 

Here Af are now the SU(2) gauge fields and g is the weak coupling constant. The 
change of baryon number AB is related to AiVcs by 

AB = n f AN cs (47) 

where nf is the number of fermion families. A change AiVcs ~ 1 requires the formation 
of a magnetic field configuration with energy of order (g 2 AR)~ l where AR is its spatial 
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extent. In order for this configuration not to be Boltzmann suppressed, its energy 
should not be larger than the temperature which requires AR >(5 ,2 T)~ 1 . On the other 
hand, the size of the field configuration cannot exceed the correlation length which is 
of order (g 2 T) -1 . Thus the length scale relevant to the problem is just the soft scale 
(g 2 T)~ 1 at which finite temperature perturbation theory breaks down. The amplitude 
of this field configuration is A ~ gT. As we have argued above, field configuration of 
this size evolve with a frequency of order g 4 log(l/g)T corresponding to a time scale 
At ~ (gr log(l/gr)T) -1 . Thus the rate for a change of baryon number B per unit time 
and unit volume can be estimated as 

T = Kg 2 \og{l/g){g 2 T)\ (48) 

where k is a non-perturbative coefficient which does not depend on the gauge cou- 
pling g. 

What might be the use of the effective theory for the soft gauge fields derived in this 
letter for non-perturbative lattice computations of real time correlators like ([I])? The 
time scale for non-perturbative dynamics is much larger than the corresponding length 
scale. Therefore the time derivative on the l.h.s. of Eq. ( |39"D can be neglected and one 
can write 

[D J ,F^(x)] = - 1 A i (x)+C(x). (49) 

This equation should be easy to implement in a lattice calculation. Note that we have 
obtained this equation by assuming that the ultraviolet cutoff \i for the soft fields is 
as small as \og(l / g)g 2 T . Therefore the lattice cutoff should not be chosen too large. 
It would be interesting to see whether the results obtained via Eq. (ffl) depend on the 
cutoff. If they do, there might be a way to include counter-terms such that the cutoff 
dependence is canceled. 

It is conceivable that, unlike the hard thermal loops 0, BD1, the effective theory for 



the soft fields in a classical lattice gauge theory has a similar structure as in the quantum 
theory. If this turned out to be the case, one could match the coefficients in the classical 
counterpart of Eqs. Q59D with the coefficients obtained in this letter. Then a numerical 
computation in the classical lattice theory should give the correct leading order result 
for a non-perturbative correlation function. The difficulty of this approach is that sub- 
leading contributions to the effective theory are only logarithmically suppressed and 
can therefore not easily be identified in numerical data. It might be possible to improve 
the matching such that also non-logarithmic terms are included. This would require 
an extension of the calculation presented here. 

In Ref. [O an algorithm was employed which, in addition to classical gauge fields, 



contains massless particles which interact weakly with the gauge fields. This approach 
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solves the problem of implementing hard thermal loops on a lattice and it appears to 
be equivalent to the starting point for the present calculation, i.e., the hard thermal 
loop effective theory described by Eqs. With this algorithm it is again difficult 

to distinguish logarithmic from non-logarithmic contributions. Probably only a combi- 
nation of the different methods discussed above can determine the rate for electroweak 
baryon number violation. 
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